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Proposal for absolute CEP measurement using 0-to-/ self-referencing 

S. B. P. Radnor, P. Kinsler,0 and G. H. C. New 

Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2AZ, United Kingdom. 

(Dated: March 17, 2009) 

We show how to adapt a — / self-referencing technique [Q, ^ to provide a single shot absolute 
Carrier Envelope Phase (CEP) measurement by using the CEP reference provided by difference 
frequency generation (DFG) between the spectral wings of the fundamental pulse. Usually, the 
beat between the input pulse and the DFG signal then provides feedback with which to stabilize 
the CEP slip in a pulse train. However, with a simple extension we can get a single shot absolute 
CEP measurement. Success relies on having well characterized input pulses, and the use of accurate 
propagation models through the nonlinear crystal - these enable us to construct a mapping between 
the experimental measurement and the CEP of the optical pulse. 



Note: The research in this paper comprised part of S. 
B. P. Radnor's PhD dissertation and was completed 
in conjunction with the other authors. It was initially 
reported in 2007 at ECLEO in Munich. This text was 
eventually finalised by P. Kinslcr. 



I. INTRODUCTION 

The management of Carrier Envelope Phase (CEP) is 
of paramount importance in attosecond physics and fre- 
quency metrology. In the absence of intervention, the 
CEP shifts from one pulse to the next in a pulse train, 
a feature that is unacceptable in experminents involving 
few-cycle pulses. The first step is to stablize the CEP, 
the second to measure the CEP offset once stabilized, 
and the third to be able to create any desired CEP to 
order. 

Recently, a — / self-referencing technique for CEP 
stabilization was developed |^ as an alternative to the 
/ — 2/ schemes previously used for CEP stabilization. 
The scheme was based on difference frequency genera- 
tion (DFG) 1^ between the spectral wings of the input 
pulse, a process in which the overall CEP is cancelled 
out. The beat between the DFG signal and the extreme 
low-frequency wing of the pulse then provides feedback 
to the source laser to enable elimination of the CEP drift. 

Measurements of the absolute CEP have been based 
on a number of techniques including photo-ionization Jq] 
high harmonic generation Q , and plasma generation 
Some innovative methods for single-shot CEP measure- 
ment based on spectral interference have also been sug- 
gested [|[ |l^ . These methods rely on CEP dependent in- 
terference occurring between various harmonics. In the 
case of Mehendale et al. ||^ this involves interference 
between the second and third harmonics, whereas the 
work by Kakehata et al. relies on interference be- 
tween a delayed fundamental and its second harmonic. 
Both mechanisms work on the basis of a relative CEP 



dependent relationship being enforced by the nonlinear- 
ity. Though the schemes are interesting, it is not clear 
how sensitive their interference assumptions are to inten- 
sity fluctuations, propagation distance, and so on. 

In the present paper, we show how the — / technique 
^ can be extended to enable a measurement of ab- 
solute CEP to be made Provided the evolution of 
the pulse within the nonlinear crystal used for the DFG 
can be accurately mapped, we show that the CEP can 
be recovered from the detailed characteristics of the beat 
signal used as the feedback source in the original experi- 
ment. We test the robustness of this technique to phase 
and intensity variations. 

In section ||, we establish a rigorous definition of CEP 
that provides a sound basis for section |l|, where we 
discover how the absolute CEP can be recovered from 
an interferometric measurement. In section we de- 
scribe the numerical techniques needed to extract the 
CEP value from the interference record, and in section 
we test their reliability. In section 0, we show how 
this enables one to measu re th e absolute CEP of an input 
pulse, followed in section VII by our conclusions. 



'Electronic address: Dr.Paul.Kinsler@physics.org 



II. CARRIER ENVELOPE PHASE 

Developments in ultrafast optical pulses have led to the 
production of sub-cycle pulses. In these limits, robust 
definitions are needed to fully characterise the pulse, as 
common descriptors can become ambiguous or fail. Per- 
haps the best example of this ambiguity is the representa- 
tion of a pulse with a carrier and envelope, where even as 
early as 1946 it was known that carrier envelope decom- 
positions (of radar pulses) were not unique pT[ |. Brabec 
and Krausz went some way towards dealing with these 
issues by suggesting a definition for the central frequency, 
and stating that an envelope definition is only valid if it 
remains invariant under a change of phase |12f| . 

The most natural way of defining CEP would appear to 
be based on a time-domain picture in which the time in- 
terval between the peak of the pulse and the closest maxi- 
mum or minimum is measured as a fraction of the optical 
period. However, for more complex pulse shapes there is 
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no unambiguous way to determine peak of the pulse^ , or 
to determine the period when the spectral bandwidth is 
broad. Further, a fixed pulse envelope in concert with 
a varying carrier phase generates a pulse that does not 
guarantee a constant energy or satisfaction of the the 
zero- force condition (see e.g. jl^). We therefore adopt a 
purely spectral approach in which we define the absolute 
spectral CEP of a pulse using the equation 

*o(^^) = (f>o + ij{oj). (1) 

The CEP we wish to measure is 00, while the relative 
spectral phase ip{uj) is set to zero at a chosen centre fre- 
quency luq. Since 4'{wo) = 0, it follows that (pQ = '^{ujq). 
We note that ip{i-o) is a measurable quantity even for ul- 
trashort pulses - It is possible to determine the relative 
spectral phases of few-cycle pulses |l5[^. This can be 
achieved to an accuracy of 0.04 rads using the SPIDER 
technique jl^. This tells us the relative CEP of all 
frequency components of the pulse, something which we 
use in the scheme presented in this paper. 

Having considered a single pulse, we now need to con- 
sider a pulse train, the spectrum of which is a comb of 
equally spaced frequency components with the lowest fre- 
quency tooth at ujcep- The phase slip between pulses is 
A = LJcepT, where r is the time interval between sucessive 
pulses in the train. Hence the absolute spectral phase of 
the n-th pulse in a train is 

*„(lj) = 0„-fV(^), (2) 

where (f>n — (/^o + nA. 

III. SCHEME 

If we know the relative spectral phase ip{uj), and we 
can determine the absolute CEP of any one frequency, 
then we can calculate the absolute CEP of any other. 
The advantage of a DFG process [|| is that it provides 
us with just such a reference CEP. 

To get efficient DFG, a nonlinear crystal (such as 
MgO:LN) is periodically poled to phase match a par- 
ticular frequency mixing process. The poling period is 
chosen so that selected frequencies on the upper (lo^) 
and lower (u;_) wings of the pulse spectrum are phase 
matched for DFG at tOd ~ w+ — w_, as illustrated on 
fig. ^ The nonlinear interaction generates a polarization 
term whose phase is 

*p = «'„(w+) - *„(w_) = ^P{lo+) - ^p{uj-). (3) 



^ It is of course possible to invent schemes for generating a suitable 
centre position, e.g. by calculating a weighted average over its 
intensity profile (see e.g. |12[ ) - however this can perform poorly 
for pulses with satellite peaks 

^ In Kobayashi et al.'s Fig. 4(b), i/> extends from 550nm to SOOnm. 
In our scheme, will most likely need to be characterised further 
down into the low frequency wing - our ui^/uiii ratio is 2.6; for 
Kobayashi et al. the ratio was only 1.5. 



Since the unknown offsets 4>o and A have canceled out, 
can be calculated from tp{uj), and provides a CEP ref- 
erence. Unfortunately, this cannot be measured directly, 
so we have to analyse how both the incident pulse and 
the DFG propagates through the crystal, and how they 
interfere and are measured. 

First, let us consider the DFG component. The polar- 
ization term (with absolute phase 5*^) will then generate 
a DFG signal that exits the crystal with a phase shifted 
by an amount Sp, so that 

'^d = "fp + Sp. (4) 
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FIG. 1: Diagram of the DFG signal being generated, compar- 
ing the phase matched DFG component to the wing of the 
input pulse. 

We can see how this shift 6p arises in the idealised 
case where a detectable DFG signal could be generated 
from a thin layer of dispersionless x'^' medium. Here 
Sp — 7r/2, since the DFG signal is just the integral of the 
driving polarization. In general, however, Sp will be a 
complicated function of both ^'p and the pulse intensity 
and profile, since it results from propagation through a 
crystal which is both nonlinear and dispersive. 

Second, we need to consider the wing of the input pulse 
spectrum at frequency u>d, which is co-propagating with 
the DFG component discussed above. This part of the 
pulse had an initial phase \E'„(cjrf), but this changes as 
it propagates through the crystal, and on exit it has a 
phase \E'^(wd) that has been shifted by Sd, i.e. 

Ki^d) = -fnioJdJ+Sd- (5) 

Although Sd will predominantly arise from the disper- 
sion, but there will also be additional contributions, due 
to other DFG processes (i.e. those not involving 0;+, 
from SPM, four-wave mixing or other processes. How- 
ever, with careful design these can be minimized, and so 
will only add an unimportant (but nevertheless calcula- 
ble) offset to Sd- 

Lastly, we apply our understanding of both the DFG 
and pulse propagation to determine the interference be- 
tween them, and what would be measured on a photode- 
tector. As the pulse exits the crystal, the DFG compo- 
nent (with phase '^d) will interfere with the wing of the 
input pulse at that same frequency, which now has phase 
^'^(wd). Thus the photodetector sees an interference be- 
tween the DFG and the spectral wing at frequency tOd- 



2 



ABSCEP 



Dr.PaTil.Kinslcr@physics.org 



From this interference it is possible to infer the relative 
phase Si between the contributions, which is 



6i = ^^(wd)-*^ 



(6) 
(7) 



In order to optimize the visibility of this interference, 
we need to ensure that the DFG signal and the ampli- 
tude of the wing of the input pulse at that frequency are 
comparable; if either is too dominant, the CEP sensitive 
modulation of the interference will be less detectable. 

In existing — / CEP stabilization experiments |j, 
the evolving interference signal resulting from the train of 
CEP-slipping pulses produces a beat signal dependent on 
the CEP slip. This beat is then used in feedback designed 
to reduce the CEP slip to zero, thus stabilizing the CEP 
of the pulses in the train to a fixed (but unknown) value 
^^{u) = (/)o + V'(w). 

The scheme presented here works because we incorpo- 
rate additional information based on knowledge of how 
the pulse propagates through the crystal. This means that 
we can predict the phase shifts Sp and 6d, at which point 
the interference measurement (i.e. of Si) turns a simple 
CEP stabilization into a CEP measurement. 

To make it clear how the CEP measurement is con- 
structed, we now show how the various phases in the 
scheme are related. We are specifically interested in the 
CEP (j)Q of the central frequency luq of the pulse; so it is 
useful to write down how the relative phases between the 
low frequency (w^) wing of the input pulse and its centre 
are 



(8) 



Now, by substituting in the preceeding collection of phase 
relationships (eqns. (||, |^, |5, ^), we can get 

*„(wo) = '/'« = ij{Lu+) ~iP{LU^) +S -ij{ujd), (9) 

where the sum of all the nonlinearity and propagation 
phase shifts is S — Sp + Si — Sd- 

Thus if we [a] know the relative phase spectrum ip^uj) 
of the input pulse(s), [b] understand the phase shifts in- 
troduced by propagation through the crystal (Sd) and of 
the DFG (Sp), and then can [c] measure the interference 
phase difference (Si), we will know every part of the RHS 
of eqn. (^ - i.e. we know the absolute CEP 0„ of the 
incoming pulse. 

The most challenging part of the CEP measurement 
is determining the Sp and Sd contributions. Fortunately, 
we can avoid having to calculate them individually by 
calculating them all at the same time in a numerical sim- 
ulation, as discussed in the following section. 



IV. MODELING 

The modeling is a crucial part of the scheme, since it al- 
lows us to determine how differing input CEPs map onto 



the detected interference measurements. After choosing 
our crystal and evaluating its parameters, and character- 
ising the pulses in our pulse train (particularly ijji^): we 
run a set of simulations over the range of CEPs. The 
results can then be used to build a map between the in- 
put CEP and the interference signal. To do this we need 
to take the spectrum of each output pulse from a simu- 
lation, and integrate over the detector response. In our 
results, we assume the spectral response of an InGaAs- 
Hamamatsu PD as used by Fuji and others ^ ||. We 
then need to check the mapping, and ensure that we will 
get the required level of discrimination between interfer- 
ence measurements from pulses with different CEP's. 

To do the propagation part of the modeling, we solve 
Maxwell's equations using the PSSD technique for a 
chosen crystal thickness. In existing — / experiments, 
the crystal thickness is typically of the order of millime- 
tres, so diffraction is negligible. We consider MgO:LN, 
with parameters taken from Further, with a suitable 
choice of nonlinear crystal (i.e. MgO:LN, and parame- 
ters from the nonlinear x*-^-* interaction occurs 
only in the extraordinary polarization (i.e. is e -|- e — > e), 
allowing us to further reduce Maxwell's equations to 



(10) 
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E + x^^^ *E + x^^^E^ + X^^^E^ ill) 



where x^^-* contains linear dispersion, and any nonlinear 
response is assumed to be instantaneous. 

In existing experiments ||, |j, a MgO:LN crystal is 
periodically poled at 11.21^m, and is optimized for DFG 
between the wings of the fundamental: a;+ (3.04 x 10^^) — 
w_(1.885 X 10^5) = Wd(1.155 X lO^^) [rad s'^]. The peak 
pulse power was calculated to be ~ 5 x lO^^W/cm^, with 
a duration of ^ 6fs (830 nm carrier). At these pulse 
powers in MgO:LN, the relative nonlinear strengths in 
the crystal at the pulse peaks are x^^^'E — 0.082 and 
T^(3)^2 _ 0.0027. This means that the self-phase modu- 
lation (SPM) distance is Lspm = 0.26mm, implying sig- 
nificant SPM over a 2mm crystal, along with other x^^-* 
effects such as 4- wave mixing, further complicating prop- 
agation and DFG. Modelling of extreme SPM (only) and 
sensitivity to CEP has been considered by Kinsler pC| ] 
and also Genty et al.pl||. 

With these parameter values, the pulses inside the 
crystal decohere rapidly because of their high intensity 
and short duration. In combination with the relatively 
small energy content within the spectral wings, the DFG 
does not continue to grow throughout the whole crystal, 
but does so only over several coherence lengths. Never- 
theless, the DFG signal can still be made large enough for 
the photodiode (PD) to detect a beat against the wing 
of the input spectrum. 

In our modeling, we retained the bulk of these parame- 
ter values to match experiment, but adjusted the crystal 
thicknesses to 100/im. This is because thinner crystals 
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still create sufficient DFG, but also have the advantage of 
producing a cleaner CEP to interference signal mapping. 



to the input CEP, the phase difference between sub- 
sequent pulses will be zero. 



V. TESTING THE CEP RESPONSE AT cua 

In order to make our scheme work, we need to guaran- 
tee that the response of the measured interference signal 
depends on the CEP in a reliable way, and that it is suffi- 
ciently insensitive to other pulse characteristics, such as 
intensity. Obviously these will depend on the particu- 
lar design of the experiment, e.g. the size of the chosen 
nonlinear crystal, the periodic-poling length, pulse fre- 
quency, and so on. In this section, we use the parame- 
ters described in the previous section to test the stability 
of the pulse propagation and interference signal against 
CEP variation and intensity fluctuations. 



A. Response to CEP slip A 

To test whether the interference signal will behave as 
expected, we did a set of simulations. Each simulations 
started with the same parameters, except for a cumula- 
tive "shot-to-shot" CEP shp of A = tt/IO. Fig. | shows 
example input pulses. 




FIG. 2: Initial pulses with ^'„(tJo) = and 7r/2 (cosine & 
sine) . 

Let us consider an idealized case, with only DFG and a 
linearly propagated input pulse present. Here we should 
see two regimes when comparing the phase spectrum of 
output pulses generated by input pulses differing by a 
CEP slip A. 



In between these two regimes will be a transition re- 
gion where the two are comparable; this is just the regime 
in which we look for the interference between the incom- 
ing pulse (with its phase slip A, as described in point 
1 above); and the phase stabilised DFG (with no phase 
slip, as described in point 2 above). Both regions, and 
the transition region of interference between them can be 
seen on fig. |^. It is important to note that, the change 
between different pairs of simulations is small - irrespec- 
tive of the absolute CEP values chosen for the two pulses. 
The only departures are at the narrow spikes caused by 
the (expected) strong CEP sensitivity near the nodes of 
the pulse. As an aside, we could (if desired) also estimate 
the linearity of the response to the CEP slip as done by 
Kinsler [ pO[ in the extreme nonlinear regime; that work 
also suggests more systematic tests of simulation pairs to 
investigate the CEP dependence. 

Here, however, we are satisfied by the fact that fig. ^ 
not only shows the predicted DFG phase stabilized region 
at low frequencies, but also the existance of a transition 
region where strongly modulated phase sensitive interef- 
erence takes place. 
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FIG. 3: Plot of the shot-to-shot changes of phase. The right 
hand side of the figure u > uii shows the expected tt/IO phase 
change, while the uj < ui region (left) represents the phase 
stable signal. The noise present at w ~ wi is the interference 
between and ujpm- The dotted lines represent the spectral 
range of the PD (colours: green=initial, red=final and blue 
are the intermediate A(j>ui curves). 



1. The spectrum of the input pulse dominates the con- 
tribution from the DFG signal, typically this oc- 
curs closer to the centre of the pulse spectrum, i.e. 
u! > uji In this case, the CEP of the pulse domi- 
nates, so that the phase difference between subse- 
quent pulses will just be the inter-pulse CEP slip 
A. 

2. The DFG signal dominates the contribution from 
the input pulse, typically this occurs for low fre- 
quencies, i.e. Lu < LOi. Since the DFG is insensitive 



B. Response to intensity fluctuations 

Our scheme relies on nonlinear interactions, but these 
are strongly intensity dependent. This means that inten- 
sity changes between pulses in the train might change the 
interference signal in a way that masks the CEP sensi- 
tivity of the interference signal. We now test the effect of 
intensity fluctuations by fixing the CEP for a set of sim- 
ulations, whilst making shot-to-shot changes in intensity 
spanning a range of ^ ±1%. 
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The results of this simulation set are displayed on fig. 
^, which demonstrates that intensity variation has a rel- 
atively weak effect in the interference region. This means 
that the mapping between photodetector signal and abso- 
lute CEP will be insensitive to the intensity fluctuations 
in the pulse train. As a result, for our parameters, we 
can disregard the effect of intensity fluctuations when it 
comes to reconstructing the CEP of an individual pulse. 
In a more extreme nonlinear case this is not always true, 
see e.g. |2^; but here the nonlinear phase shifts (e.g. 
those due to SPM) are not strongly intensity dependent. 
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FIG. 4: Plot of for intensities ranging from 99% to 101%. 
All pulses began with (f)cE ~ 0. The red and green curves 
represents the lowest and highest intensities respectively. The 
2% range is divided into 10 equal sections (blue). 



VI. MEASURING THE ABSOLUTE CEP 

In the previous section we demonstrated that (or our 
chosen parameters) not only was the interference in the 
ujd DFG region sensitive to CEP in a controllable way, 
the effect would not be masked by intensity fluctuations. 
Since we see this clear dependence on CEP, it is possi- 
ble to determine the absolute CEP from the interference 
signal - as long as the pulse intensity and crystal param- 
eters are appropriately matched. For example, the CEP 
dependence is better behaved at some distances than at 
others - for our chosen parameters, it happens that dis- 
tances of ~ 50, lOO^m give good results. 

Fig. H shows the CEP dependent structure at a propa- 
gation distance of 50/im, where each curve is the intensity 
for a different input CEP. Using this, we can then inte- 
grate that spectral behaviour over the response of the 
photodetector to generate our mapping. The mapping 
corresponding to fig. ^ and our chosen photodiode is 
shown on fig. ^. 

However, because each interference signal value is not 
unique, we have only determined the CEP to within tt. 
To complete the determination of the CEP to within a 
2tt range we need to take two such measurements under 
slightly different conditions, e.g. using different propaga- 
tion lengths. 

Summarizing, for this CEP measurement scheme to 
succeed, we must have accurate knowledge of the follow- 
ing three things: 
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FIG. 5: Plot of the spectral intensity at uj : 



for a prop- 



agation distance of 50/im. The CEP dependent structure can 
clearly be seen, where the lowest curve (yellow) corresponds 
to 4>CE ~ and the upper (turquoise) curve corresponds to 
4'CE = 97r/10. The range is divided into intervals of vr/lO 
rads. 
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FIG. 6; Mapping between photodetector signal and pulse 
CEP: By integrating the power falling on the PD within the 
detectable range, a CEP varying signal was constructed. In 
this case = maps to the first point on the PD signal (an 
extra point has been plotted for (pcE ~ tt). 



1. The phase spectrum of the pulse. In 

eqn.(^ we see that it depends on the sums or dif- 
ference of four values of ipijjj), thus compounding 
any uncertainty in its determination. 

2. The intensity fluctuations of the pulse. The 

intensity must be controlled sufficiently well so that 
the SPM or XPM effects are minimized, and do 
not significantly alter the propagation of the pulse 
through the crystal. 

3. Pulse propagation. This can be done easily with 
PSSD simulation code (or similar), but we need 
accurate information on the initial conditions of the 
pulses in the pulse train. 

Remarkably, this can be done with a small extension 
to current — / self-referencing methods, which are cur- 
rently only used to stabilize the CEP. The extension is to 
numerically model the propagation of the pulse through 
the nonlinear crystal in order to determine the mapping 
between the photdetector signal and the input CEP. 
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VII. CONCLUSION 

We have demonstrated how an ordinary — / self- 
referencing scheme can be easily extended to measure 
absolute CEP, rather than just being used as a CEP sta- 
bilization tool. The scheme relies on a phase stable sig- 



nal being passively produced through DFG, and does not 
require strong field physics to operate. Instead, we pro- 
pose to numerically model the propagation of the pulse 
through the nonlinear crystal; and to use the information 
gained to determine the mapping between the detected 
interference signal and absolute CEP of the input pulse. 
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